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The lattice QCD simulations of charmonia and exotic charmonium-likc 
states are reviewed. I report on the first exploratory simulation of charmo- 
nium resonances above open charm threshold which takes into account the 
strong decay. The puzzles related to the first-excited scalar charmonia are 
discussed. Evidence for A"(3872) is presented along with investigation of 
its Fock components. The Zf (3900) seems to emerge from the HALQCD 
approach as a result of the coupled channel effect J/fm — DD . The 
indication for the pentaquark bound state rj c N is presented. 
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1 Introduction 


Experiments discovered a number of interesting charmonium-like states which can 
not be reconciled with conventional cc. The prominent examples are tetraquarks Z p 
with flavour content ccdu and pentquarks P c with ccuud [lj. There is no doubt that 
establishing the conventional and the exotic charmonia from ab-initio lattice QCD 
simulations is desired. 1 review state of art calculations along these lines. 


2 Lattice methodology 

The physics information on a hadron (below, near or above threshold) is commonly 
extracted from the discrete energy spectrum in lattice QCD. The physical system 
for given quantum numbers is created from the vacuum Q) using interpolator Oj 
at time t = 0 and the system propagates for time t before being annihilated by 
0{. To study a charmonium or a charmonium-like state with given J p one can, for 
example, use (9 ~ cTc, two-meson interpolators O = {cTiq)(qT 2 c), (cTic)(qr 2 q) and 
O = [cr 1 g][cT 2 g] with desired quantum numbers. After the spectral decomposition 
the correlators are expressed in terms of the energies E n of eigenstates | n) and their 
overlaps Z 7 j 


Cijit) = {n\Oi(t)o](o)\n) = Y W e ~ Ent > = (n\Oi\n ). (1) 

n 

The correlators are evaluated on the lattice and their time-dependence allows to 
extract E n and Z l n M. 

The energy eigenstates | n) are predominantly ”one-meson” states (e.g. x c i( 1P)) 
or predominantly ’’two-meson” states (e.g. DD ) - in interacting theory they are 
mixtures of those. Two-meson states have a discrete spectrum due to the periodic 
boundary condition on the finite lattices. If they do not interact, then the mo¬ 
menta of each meson is p — ^N with N 6 N 3 , and the non-interacting energies of 
Mi(p)M 2 (— p) are E n - 1 - = E\(p) + E 2 (jp) with Ei t2 {p) = (mf j2 + p 2 ) 1 ^ 2 - The energies 
E n extracted from the lattice are slightly shifted in presence of the interaction and the 
shift provides rigorous information on the scattering matrix, as discussed in Section 
[5j In experiment, two-meson states correspond to the two-meson decay products with 
a continuous energy spectrum. 

The Wick contractions for charmonium (like) states contain terms where charm 
quark propagates from source to sink, and terms where charm quark annihilates. 
The charm quark annihilation is omitted in all lattice simulations reported in this 
talk. One expects very small influence from charm annihilation on the energies of 
eigenstates of interest, but this needs to be verified in the future. It will be very 
challenging to go beyond this approximation due to a large number of light single 


1 


and multi-hadron states with the same quantum numbers. For some channels (for 
example DD scattering with 1 = 0), there are also Wick contractions where u,d,s 
annihilate and these Wick contractions have been consistently taken into account in 
all the described simulations. 



Charmonium low spectrum 



Figure 1: Spectra of charmonia below open charm threshold [Tj (left) and [5] (right). 


3 Charmonia well below open charm threshold 

The masses rrii = Ei(P = 0) are extracted from the energies obtained with cc inter¬ 
polating fields, which are extrapolated a —» 0, V —>■ oo and rn q —>■ m yhys . Since charm 
annihilation is omitted, there is no complication arising in relation to the multi-hadron 
states. Some of the recent spectra m 12 la E| are shown in Fig. [I] The results are in 
impressive agreement with experimental masses and the main remaining uncertainty 
is due to the omission of charm annihilation. 


4 Excited charmonia within single-meson approach 

The most extensive spectra of the excited charmonia have been calculated within the 
so-called single-meson approach Hadron Spectrum Collaboration in 2012 [8]. The 
continuum J PC was reliably identified using the advanced spin-identification method. 
Multiplets of hybrid states were also found and some of them carry exotic J PC . The 
single-meson treatment ignores strong decays of resonances and threshold effects. It 
gives valuable reference spectra, but can not give reliable conclusions on the near¬ 
threshold exotic states, for example. 
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Figure 2: Left: The spectrum of the vector charmonia from [3]: the diamond denotes 
the resonance mass of 0(3770), while the triangle denotes the pole mass of the bound 
state '0(25); both are obtained from DD scattering matrix. Right: The location of 
X(3872) with 1 = 0 which emerges as shallow bound state in DD scattering at 
m w ~ 266 MeV [10] and m n ~ 310 MeV [II j (update of [I2j). 


5 Scattering approach to resonances and bound 
states 


In the energy region near or above threshold, the masses of bound-states and res¬ 
onances have to be inferred from the infinite-volume scattering matrix of the one- 
channel (elastic) or multiple-channel (inelastic) scattering. The simplest example 
is a one-channel elastic scattering in partial wave l , where the scattering matrix is 
parametrized in terms of the phase shift 5i(p) and satisfies unitarity SS^ = 1 

S(E) = e 2ii ‘< E > , S(E) = 1 + 2 iT(E) , T(E) = ! . . (2) 

cot {di(E)) -1 

Liischer has shown that the energy E of two-meson eigenstate in finite volume L 
gives the elastic phase shift 5(E) at that energy in infinite volume [13]. This relation 
and its generalizations are at the core of extracting rigorous information about the 
scattering from the recent lattice simulation. It leads 5(E) only for specific values of 
E since since spectrum of two-meson eigenstates is discrete. The 5(E) or T(E ) ([2]) 
provide the masses of resonances and bound states: 


• In the vicinity of a hadronic resonance with mass rriji and width T, the cross 
section cr oc \T(p)\ 2 has a Breit-Wigner-type shape with 5(s = m 2 R ) = | 


T(p) 


-y/i r(p) 

s — m 2 R + iy/s T(p) 


1 

cot 5(p) — i 


r(p) = g 2 - 


21+1 


S 


(3) 
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The fit of 5i(p) renders m B and g or T. It is customary to compare g rather 
than T to experiment, since T depends on the phase space. Lusher’s approach 
has been verified on several conventional mesonic resonances. 

• The bound state (B) in Mi M 2 scattering is realized when T(p) has a pole at 
Pb < 0 or p B = i\p B \ 

T = 1 - 7 = 00 , cot (d(p B )) = i, m B = E Hi (pb) ± E H2 (p B ) • (4) 

cot {Oi(p B ))-1 

The location of an s-wave shallow bound state can be obtained by parametrizing 
So near threshold and finding p B which satisfies cot (S(p B )) = i. 


6 Vector and scalar resonances 

Until recently, all charmonia above open-charm threshold were treated ignoring the 
strong decay to a pair of charmed mesons. The first exploratory simulation aimed 
at determining the masses as well as the decay widths of these resonances was pre¬ 
sented in [2j. The Breit-Wigner-type fit of the DD phase shift in p -wave leads to the 
resonance mass and width of 0(3770) (see Fig. [6] and Table [TJ. The 0(2S') in Fig. [6] 
appears as a bound state pole below threshold. 


0(3770) 

m res [MeV] 

g (no unit) 

Lat 

(m 7 r = 266 MeV) 

3774 ± 6 ± 10 

19.7 ± 1.4 

Lat 

(m n = 156 MeV) 

3789 ± 68 ± 10 

28 ±21 


Exp. 

3773.15 ±0.33 

18.7 ± 1.4 


Table 1: Parameters of the resonance "0(3770) in DD scattering. The width is 
parametrized in terms of the coupling g as T = g 2 p 3 /(6ns) [9]. 

I 11 the scalar channel, only the ground state y c o(lP) is understood and there is no 
commonly accepted candidate for its first excitation y c0 (2 P) although PDG assigns 
it to A(3915). The analytic studies loj argue that the A(3915) can probably not 
be identified with the Xco( 2P), while [IB] finds that the experimental data on A"(3915) 
could actually be compatible with J pc = 2 ++ . These studies suggest that a broad 
structure observed in the DD invariant mass represents y c o(2 P). The phase shift for 
DD scattering in s-wave was extracted from lattice [9], but the resulting phase shift 
also does not allow a clear answer to the puzzles in this channel. The lattice data 
provide an indication for a yet-unobserved narrow resonance slightly below 4 GeV 
with r [y( 0 —>■ DD] below 100 MeV. A scenario with this narrow resonance and a pole 
in the DD scattering matrix at y c o(lP) agrees with the energy-dependence of the 
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phase shift [9], Three other scenarios are not supported by the lattice data: just one 
narrow resonance, just one broad resonance (proposed in Guo & Meissner [FT] and 
Olsen [E]), or one narrow and one broad resonance. Further experimental and lattice 
QCD efforts are required to map out the s-wave DD scattering in more detail. 
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Figure 3: Left: Energies of eigenstates in Z p channel from simulation [T7| are given 
by black circles. They correspond to the expected two-meson states, whose non¬ 
interacting energies are given by the horizontal lines. Right: Binding energies of ?/ c 
and nuclei with A nucleons from [115]. The bound state at A = 1 is a pentaquark 
candidate 7j c N. 


7 Charmonium-like states X(3872) and F(4140) 

The A"(3872) lies experimentally on D°D threshold and its existence on the lattice 
can not be established without taking into account the effect of this threshold. This 
was first done by simulating DD scattering in [TO] , where a pole in DD scattering 
matrix was found just below the threshold in I(J PC ) = 0(1 ++ ) channel ([!]). The 
pole is associated with a bound state X(3872) and its location is shown in Fig. [6} 
The more recent simulation using HISQ action confirms the existence of the pole just 
below the threshold pa nu. 

The lattice study [T9] investigated which Fock components are essential for ap¬ 
pearance of A"(3872) with / = 0 on the lattice. The energy eigenstate related to 
A(3872) appears in the simulation only if DD as well as cc interpolating fields 
are employed. The X(3872) does not appear in absence of cc interpolators, even if 
(localized) interpolators \cq\^ c [cq\^ c or [cg]6 c [cg]g c are in the interpolator basis. This 
indicates that cc Fock component is most likely more essential for A"(3872) than the 
diquark-antidiquark one. 
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The search for charged X (3872) was performed in J pc = 1 ++ channel with DD *, 
J/ipp, [cd] 3c [cu)^ c and [cd] 6c [cu]g c interpolators and no candidate was found [19]. The 
reliable search for the neutral 1=1 state would need to incorporate isospin breaking 
effects [20], but that has not been performed on the lattice yet. 

The experimental candidate Y( 4140) with hidden strangeness was observed in 
J/ip(p invariant mass and has unknown J p at present. The lattice search was per¬ 
formed in J pc = 1 ++ channel with D S D S1 J/ip(p, [cs][cs] interpolators and no candi¬ 
date was found |IH|- The s-wave and p-wave J/ ip <p scattering phase shift from [2TJ 
do not support the resonant structure. 
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Figure 4: The HALQCD results related to Z+(3900) where scattering amplitude / 
is defined through dcr/dfl = \ f(E )\ 2 [22j. Top left: potentials V(r) for three channels. 
Top middle: crossed-channel potentials V(r) with sizeable J/ipn — DD . Top right: 
Poles of the scattering matrix in the complex energy plane. The quantity related to 
the number of particles Irn [/,;] oc for channels i = J/ipn, DD compared to the 
experimental shapes for J/ipn [23] and DD [24j. 
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8 Charged charmonium-like states Z c 

The lattice search for the manifestly exotic states Z F with flavour content ccdu and 

jG(jPC) = !+(!+- 

') is very challenging since the experimental candidates lie above 
several thresholds and can decay in several final states via strong interaction. The 
reliable treatment requires simulation of coupled channels and extraction of coupled 
channel scattering matrix. This has been done using the Liischers method only for 
one system in the light sector [25]. 

The lattice search for Z F in the energy region below 4.3 GeV m has therefore 
been performed using a simplified approach. The challenge is that there are 13 two- 
meson eigenstates (J/07T, DD , rj c p, ...) with I G (J FC ) = 1 + (1 H ) in this energy 
region for the L ~ 2 fm, as shown by horizontal lines in Fig. [3} The explicit 
simulation of all these coupled channels (with a large number of meson-meson and 
diquark-antidiquark interpolators) renders 13 expected two-meson eigenstates in the 
region (shown by black circles). No additional energy eigenstate has been found, 
therefore no candidate for Z F . This has been recently confirmed also by a simulation 
using H1SQ action [12]. Note that an extra eigenstate (in addition to expected two- 
meson states) has been found for all resonances and all bound states that have been 
established on the lattice up to now (this implies, for example, also for the bound state 
A"(3872) and resonance -0(3770) discussed in this talk). The absence of an additional 
eigenstate for Z c (3900) could indicate that this experimental state is not related to a 
conventional resonance pole, but could be a manifestation of coupled channel effect. 

This possibility was investigated within the so-called HALQCD approach [26] to 
extract the scattering matrix of the coupled channels DD , J/ipir and p c p [22]. The 
potential Krj/p— ti-j/pCQ related to Nambu-Bethe-Salpeter equation is determined be¬ 
tween the J/0 and 7r as a function of their separation r. This is presented in Fig. [4] 
together with potential for the other two channels, as well as potential corresponding 
to coupling between different channels according to coupled-channel HALQCD for¬ 
malism [27]. The off-diagonal potential between channels 7 tJ/i/j and DD is larger 
than other potentials, which seems to indicate a sizeable coupled channel effect near 
DD threshold. The 3x3 matrix of potentials renders 3x3 scattering matrix[27], 
This in turn gives a quantity related to the number of events iVj oc cr,j,; in each 
channel i as a function of invariant mass, which is compared to experiment in Figure 
4 Both and N DD * show enhanced peaks above DD threshold that resemble 

experimental line shapes for Z+( 3900). The pole structure that emerges from the 
extracted scattering matrix is shown in the complex energy plane in Fig. [4] Two 
poles are found at Im[E cm \ < 0 and Re[E cm ] < run + m^*, while a pole would be 
expected at Re[E cm ] > m D + if Z c (3900) would be conventional resonance. The 
authors conclude that Z F (3900) most likely does not correspond to the conventional 
resonant pole, while the cross section peak is related to sizeable 7r J/-0 — DD coupled 
channel effect, ft remains to be explored whether the extracted scattering matrix is 
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in agreement with the absence of an additional energy eigenstate near 3.9 GeV in 
mm- The conclusions based on the HALQCD approach need to be verified using 
the Liischer-type approach, which has been shown to render reliable rigorous results 
for several conventional meson resonances. The HALQCD approach has not been 
thoroughly tested on the conventional meson resonances yet. 

The s-wave and p-wave phase shifts near DD threshold were determined using 
only DD interpolating fields in [28], which may not be reliable since the ground state 
of the system is The authors conclude that no evidence for Z+(3900) is found. 


9 Pentaquarks and conucleus bound states 

The NPLQCD collaboration finds an interesting evidence for a r] c N bound state 
approximately 20 MeV below rj c N threshold, as shown in Fig. [3] [18] (N denotes 
nucleon). To my knowledge, this is the only pentaquark candidate containing cc from 
lattice studies up to now. As the simulation is done at 577(3) flavour symmetric point 
corresponding to m n ~ 800 MeV, it is not clear yet whether this bound state persists 
to physical m n . LHCb has recently found two pentaquark resonances in J/ijjp, located 
about 400 MeV above threshold fT]. The lattice simulation of those is much more 
challenging due to several open channels and has not been performed yet. 

Bound states of r] c and nuclei ware also found in [18], where binding energies are 
given in Fig. [3] 
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Figure 5: The central and spin dependent potential V(r) = Vc(r ) + Vg(r)si • s 2 
between c and c for various temperatures [29 ]. 


10 Charmonium potential at nonzero temperature 

The potential V(r) between c and c was extracted using HALQCD method [2J3J at 
a range of finite temperatures T [29]. The central and spin dependent parts V(r) = 
V c (r) + Vs(r)si • s 2 are presented in Fig. [5j The lines in Vc present Cornell and 
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Debye screening fit of the results, and the resulting Debye screening mass acquires 
a nonzero value roughly at T c . The spin-dependent part has a repulsive core which 
roughly resembles the familiar parametrisation with the 8(r) function. 


11 Summary 

Recent lattice QCD studies of charmonium and charmonium-like states were reviewed. 
The main challenge for the future lattice simulations is the extraction of the scatter¬ 
ing matrix relevant to the experimentally interesting states. This will, for example, 
involve two or more coupled channels for tetraquarks Z+ or pentaquarks P c . 
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